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The spectral function A(q, uS) of doped t — J ladders is presented on clusters with up to 2 x 20 
sites at zero temperature applying a recently developed technique that uses up to ~ 6 x 10 6 rung- 
basis states. Similarities with photoemission results for the 2D cuprates are observed, such as the 
existence of a gap at (ir, 0) near half-filling (caused by hole pair formation) and flat bands in its 
vicinity. These features should be observable in ARPES experiments on ladders. The main result 
of the paper is the nontrivial evolution of the spectral function from a narrow band at x = 0, to a 
quasi-noninteracting band at x > 0.5. It was also observed that the low-energy peaks of a cluster 
spectra acquire finite line-widths as their energies move away from the chemical potential. 

PACS numbers: 74.20.-z, 74.20.Mn, 75.25.Dw 



Copper-oxide ladder compounds are currently under 
much investigation [Q. Among their interesting proper- 
ties are a spin-liquid ground state in the undoped limit, 
and the existence of superconductivity upon hole dop- 
ing ^|J|] . Recently, the first angle-resolved photoemission 
(ARPES) studies of ladder materials have been reported. 
Both the doped and undoped ladder Sri4Cu2404i 
have been analyzed, finding one-dimensional metallic 
characteristics [Q. Studies of the ladder compound 
Lai_ x Sr x Cu02.5 found similarities with La2- x Sr x Cu04, 
including a Fermi edge Core-level photoemission 

experiments for (La, Sr, Ca)i4Cu2404i documented its 
chemical shift against hole concentration Note that 
the importance of ARPES studies for other materials 
such as the high-Tc cuprates is by now clearly estab- 
lished . Using this technique the evolution with doping 
of the Fermi surface has been discussed J8| , including the 
existence of flat bands near momenta (0, w) — (tt, 0) M. 

These plethora of experimental results for the cuprates 
should be compared against theoretical predictions. 
However, the calculation of the ARPES response even 
for simple models is a formidable task. The most reliable 
computational tools for these calculations are the Exact 
Diagonalization (ED) method, restricted to small clus- 
ters, and the Quantum Monte Carlo (QMC) technique 
supplemented by Maximum Entropy, limited in doped 
systems to high temperatures due to the sign problem. 
Currently, on ladders dynamical properties can be ex- 
actly calculated at all densities only on clusters of size 
2x8 while the QMC technique in the realis- 

tic regime of large U jt (Hubbard model) has been ap- 
plied on 2 x 16 lattices only at half-filling and with 1 
hole fLU , the latter using an anisotropic ladder since for 
the isotropic case the sign-problem is severe. 

Due to the limitations of these techniques an impor- 
tant issue still unclear is the evolution of the one-particle 
spectral function between the undoped limit, dominated 
by antiferromagnetic (AF) fluctuations both on ladders 
and planes, and the high hole-density regime where those 
fluctuations are negligible. While both extreme cases are 



properly treated by previously available numerical meth- 
ods, the transition from one to the other as the hole den- 
sity x grows is still unknown. This evolution is expected 
to be highly nontrivial. For instance, the presence of 
hole-pairs in lightly doped ladders suggests the opening 
of a gap in ARPES, similar to the pseudogap of under- 
doped high-Tc cuprates |l5). Shadow-band features in 
undoped ladders jnj, which are absent at higher hole 
densities, adds to the complexity of this evolution. 

Motivated by this challenging problem, in this paper 
the density evolution of the spectral function A(q,ui) of 
doped 2-leg t — J ladders is presented. The calculation 
is carried out at zero temperature on clusters with up to 
2 x 20 sites, increasing by a substantial factor the current 
resolution of the ED techniques. These intermediate size 
clusters were reached by working with a small fraction of 
the total Hilbert space of the system jl6| . The method 
is variational, although accurate as shown below. The 
improvement over previous efforts lies in the procedure 
used to select the basis states of the problem [|l7|. The 
generation of the new basis is in the same spirit as any 
technique of the renormalization-group (RG) family. If 
the standard 5 2 -basis is used (3 states per site), expe- 
rience shows that a large number of states is needed to 
reproduce qualitatively the spin-liquid characteristics of 
the undoped ladders. The reason is that in the 5 z -basis 
one of the states with the highest weight in the ground 
state is still the Neel state, in spite of the existence of a 
short AF correlation length £,af- A small basis built up 
around the Neel state incorrectly favors long-range spin 
order. However, if the Hamiltonian of the problem is ex- 
actly rewritten in, e.g., the rung-basis (9 states/rung for 
the t—J model) before the expansion of the Hilbert space 
is performed, then the tendency to favor a small £af is 
natural since one of the dominant states in this basis 
for the undoped case corresponds to the direct product 
of singlets in each rung, \S), which has £af = along 
the chains. Fluctuations of the Resonant- Valence-Bond 
(RVB) variety around \S) appear naturally in this new 
representation of the Hamiltonian leading to a finite £af- 
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Note that \S) is just one state of the rung-basis, while 
in the S^-basis it is represented by 2 Nr states with N r 
the number of rungs of the 2-leg ladder. In general a few 
states in the rung-basis are equivalent to a large number 
of states in the .S^-basis. Expanding the Hilbert space in 
the new representation is equivalent to working in the S z - 
basis with a number of states larger than can be reached 
directly with present day computers. Here for simplic- 
ity this technique will be referred to as the Optimized 
Reduced Basis Approximation (ORBA) p7||. 




FIG. 1. ecs vs x for the t 



4 5 
d 

J model using the method 



described in this paper for the two clusters indicated (t = 1). 
DMRG results with OBC are also provided for comparison; 
(b) Ground state hole-hole density correlations vs distance for 
a variety of clusters and densities as indicated. Some results 
using DMRG with PBC are also shown; (c) Ground state 
staggered spin-spin correlations vs distance d (= |i — j|) along 
the leg opposite to where site i is located. Number of holes are 
indicated. Some results with DMRG are also shown. Open 
(full) symbols are for 2 x 16 (2 x 20) clusters. 

As a first step, let us compare ORBA predictions for 
equal-time observables against DMRG results (l8) for the 
same clusters. Here a coupling J/t = 0.4 is used. Its par- 
ticular value is important: if J/t is smaller, then pairs 
are lost while if it is larger superconducting correlations 
are important. Only in a small window of J/t is that 
the ground state can be considered as formed by weakly 
interacting hole pairs, a regime that we want to inves- 
tigate in this paper for its possible connection with the 
phenomenology of high-Tc at finite temperature. Fig. la 
contains the ground state energy per site ecs y s x using 
~ 2 - 3 x 10 6 states in the rung-basis. The DMRG ener- 
gies are obtained with m = 200 states and open boundary 
conditions (OBC). Both sets of data are in good agree- 
ment On 2 x 20 clusters, the rung-basis approach al- 
lowed us to study up to 6 holes |2(]] which has a full space 
of ~ 10 14 states (for zero momentum and total spin), 



while the largest previously reported exact study on a 
2 x 10 cluster and 2 holes needs a~5x 10 5 basis fH)|,|llJ. 
Calculating the binding energy, or the chemical poten- 
tial ji vs x, from ecs supplemented by the energies for 
an odd number of electrons, a tendency to pair forma- 
tion at low hole-density was observed [ |10|]2l|| . Fig. lb 
contains the hole-hole correlations at several densities, 
compared (in one case) with PBC DMRG results. In 
Fig.lc spin-spin correlations are shown. The rung-basis 
properly reproduces the existence of a small £af in the 
ground state, that decreases as x grows. Size effects are 
not large for the clusters studied here, and good agree- 
ment with DMRG is observed. This technique captures 
the essence of the ground state behavior p2[. 
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FIG. 2. A(q,u) for the bonding band q y = on a 2 x 20 
cluster, (a) corresponds to x = and (b) to x — 0.10. The 
insets correspond to 2 x 16 clusters. 

To produce dynamical results Ref. |17|] was followed, 
namely ~ 10 — 20% states of the reduced basis TV-holes 
ground state \ipo) were considered [^3| and the reduced 
subspace with, e.g., N + 1 holes and q-momentum was 

obtained through d q |Vo) (0\ = Ej eiqj Cj' with 4 the 
hole creation operator at j, dropping the spin-index, in 
the rung-basis). All states generated by this procedure 
were kept, and we worked in such a subspace in the subse- 
quent iterations of the continued fraction expansion J24] 
Only the bonding band subspace is discussed here p5[ . 
The (^-functions have a width O.lt throughout the paper. 

Fig. 2a corresponds to the undopcd limit. A sharp peak 
is observed at the top of the PES spectra, maximized at 
momenta q x — 77r/10, i.e. close to the Fermi momen- 
tum for noninteracting electrons — 0.667T. The band 
defined by those peaks has a small bandwidth, as in 2D 
models, due to the interaction of the injected holes with 
the spin background |24|,|2(|. Note that all peaks at mo- 
menta q x > 7r/2 carry a similar weight and the disper- 
sion is almost negligible. This unusual result is caused by 
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strong correlation effects. The PES weight above q£ , e.g. 
at q x = tt, is induced by the finite but robust £af, and 
its existence resembles the antiferromagnetically induced 
"shadow" features discussed before in 2D models j27j . 

Fig. 2b contains results at low but finite hole-density. 
Several interesting details are observed: (i) the PES band 
near /x is flat. This should be an ARPES observable re- 
sult resembling experiments in 2D cuprates, and it adds 
to the growing evidence linking the physics of ladders and 
planes; (ii) q x — tt (ir/2) PES has lost (gained) weight 
compared with x = 0; (iii) the total PES bandwidth has 
increased; and (iv) the IPES band is intense near q x = tt, 
and it is separated from the PES band by a gap. The 
observed gap is A ~ OAt and it is caused by hole pair- 
ing. The DMRG/PBC binding energy calculated for the 
same cluster and density is ~ 0.32i [m — 200, trunca- 
tion error ~ 10~ 4 ). In the overall energy scale of the 
ARPES spectra, this difference is small and does not af- 
fect the study of the evolution of the dispersion shown 
here. Note that the results of Fig. 2b are similar to those 
observed near (tt, 0) using ARPES in the 2D cuprates 
normal state [flO . 
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FIG. 3. Same as Fig. 2 but using a 2 x 16 cluster and the 
hole densities and momenta indicated. In (a) PES, the results 
are also shown using a width 0.001£ (and a different vertical 
scale) to visualize the individual <5-functions. The number 
next to each broad peak is the number of poles contributing 
to it (some are difficult to resolve due to their small weight). 
In (b) and q = (n/2,0) results obtained using 1.1, 2.8, 4.5 
million states are shown (dashed lines) from left to right, to 
illustrate the convergence. The solid lines were obtained with 
6.0 million states. At this density the Hilbert space is maxi- 
mized for the 2 x 16 cluster. 

Fig. 3a contains results at x = 0.1875. The trends ob- 
served at x — 0.1 continue, the more dramatic being the 
reduction of the q x — tt PES weight caused by the de- 
crease in £af ■ The lost weight appears in the q x = tt 



IPES signal. The gap is still observed in the spectrum. 
Weak BCS-like features both in PES and IPES near q£ 
can be seen. Fig. 3b contains data at x = 0.3125, and 
up to 6 x 10 6 states. The Hilbert space is maximized at 
this density for the 2 x 16 cluster. Now the result resem- 
bles more a noninteracting system on a discrete lattice. 
The IPES signal is no longer very flat, and the IPES 
band now has a clear energy minimum near the momen- 
tum where PES is maximized. Fig. 3c contains results 
for x = 0.5 where a quasi-non-interacting dispersion is 
obtained using about 3 x 10 6 states in \ipo). The inset 
shows that the trend continues at lower electronic densi- 
ties. The bandwidth evolves from being dominated by J 
near half-filling, to having t as natural scale at x ~ 0.3 or 
larger. This evolution is smooth, yet nontrivial, following 
the reduction of ^af with doping. 

A conceptually interesting issue in the context of finite- 
cluster spectra of electronic models is whether finite line- 
widths for the dominant peaks can be obtained by such 
a procedure. Studying the small clusters reached by ED 
techniques it naively seems that those peaks are usually 
generated by just one (^-function (one pole). However, in 
the bulk limit, peaks away from the Fermi level should 
have an intrinsic width. How can we reach such a limit 
from finite clusters? One possibility is that as the cluster 
grows, the number of poles N p in a small energy win- 
dow centered at the expected peak position must grow 
also, with their individual intensities becoming smaller 
such that the combined strength remains approximately 
constant. While this idea seems reasonable, it still has 
no explicit verification, but the intermediate size clusters 
reached in this study allow us to test it. Consider as 
an example Fig. 3a where the actual energy and intensity 
of the poles contributing to the main features are shown 
explicitly. As the peaks move away from the top of the 
PES band, N p was indeed found to increase providing ev- 
idence compatible with the conjecture made above P8| . 

Fig. 4a contains the main-peak weights in the PES band 
vs density. Size effects are small. The weight at q x = tt 
diminishes rapidly with x, following the strength of the 
spin correlations of Fig.lc. Overall the region affected 
the most by spin correlations is approximately x < 0.25. 
Fig. 4b summarizes the main result of the paper, pro- 
viding to the reader the evolution with x of the lad- 
der dominant peaks in A(q, u>). The area of the circles 
are proportional to the peak intensities. At small x a 
hole-pairing-induced gap centered at fx is present in the 
spectrum, both the PES and IPES spectra are flat near 
(tt, 0), and the band is narrow. The PES flat regions at 
high momenta exist also in the undoped limit, where they 
are caused by the short-range spin correlations. Actually 
the resolution in densities and momenta achieved in this 
study allow us to reach the conclusion that the undoped 
and lightly doped regimes are smoothly connected. As 
x grows to ~ 0.3, the flat regions rapidly loose intensity 
near (7r,0), and the gap collapses. 
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Summarizing, the bonding-band spectral function of 
the 2-leg t — J model has been calculated, and results 
can be used to guide future ARPES experiments for lad- 
der compounds. These experiments should observe flat 
bands and gap features near (tt, 0) in the normal state. 
The data was found to be remarkably similar to experi- 
mental results for the 2D cuprates along the (0, 0) — (it, 0) 
line. A common explanation for these features was pro- 
posed. Finally, note that the ORBA method discussed 
here introduces a new way to calculate dynamical proper- 
ties of spin and hole models on intermediate size clusters. 
The method can be applied to a variety of strongly cor- 
related electronic models. 
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FIG. 4. (a) Weight of the low-energy dominant feature in 
A(q, a)) for two momenta q x in the bonding band vs x. Cluster 
sizes are indicated; (b) Evolution with doping of the dominant 
feature band (t — 1). The open (full) circles are centered at 
the peak energies e qp below (above) fi. The hole densities are 
indicated. The area of the dots is proportional to the weight 
of the peak. The results can apply to 2D systems along the 
(tt, 0) - (0, 0) line. 

The many similarities between ladders and planes dis- 
cussed in previous literature f| suggest that our results 
may also be of relevance for 2D systems along the line 
(0,0) — (?r,0). For instance, the abnormally flat regions 
near (tt, 0) (Fig. 2b) are similar to ARPES experiments 
data for the 2D cuprates ||, and they should appear in 
high resolution photoemission experiments for ladders as 
well. Note that in the regime studied with pairs in the 
ground state, the flat bands do not cross fi with doping 
but they simply melt. When x is between 0.3 and 0.4, a 
quasi- free dispersion is recovered. The results of Fig. 4b 
resemble a Fermi level crossing at x ~ 0.3 and beyond, 
while at small hole density no crossing is observed. It 
is remarkable that these same qualitative behavior ap- 
peared in the ARPES results observed recently in under- 
doped and overdoped LSCO |2S|] . These common trends 
on ladders and planes suggest that the large energy scale 
(LES) pseudogap (~ Q.2eV) of the latter Q may be 
caused by similar long-lived d-wave-like tight hole pairs 
in the normal state as it happens in doped ladders, where 
the pairing is caused by the spin-liquid RVB character of 
the ground state || . In the 2D case a similar effect may 
originate in the finite ^af observed in the undcrdoped 
finite temperature regime (although there is no clear evi- 
dence in the 2D layered materials of a spin-liquid ground 
state). A consequence of this idea is that the pairs and 
LES ARPES pseudogap are correlated and they should 
exist as long as £af is non- negligible |3l]], a prediction 
that can be tested experimentally. 
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